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HEAT TRANSFER IN THERMAL ENTRANCE REGION WITH LAMINAR SLIP FLOW 
BETWEEN PARALLEL PLATES AT UNEQUAL TEMPERATURES 
by Robert M. Inman 
Lewis Research Center 

SUMMARY 

An analysis was made of the forced- convection heat-transfer characteristics for 
fully developed incompressible laminar flow between parallel walls with unequal tempera- 
tures under the condition that the gas density is low enough to permit a velocity slip and a 
temperature jump at the walls. The results apply along the entire length of the channel. 
The solutions contain series expansions and analytical expressions for the complete sets 
of eigenvalues and eigenfunctions. The results give the wall heat-flux requirements and 
Nusselt numbers for various values of the rarefaction parameters and dimensionless 
entrance-temperature parameter. The results show that, in general, the slip-flow heat- 
flux requirements and Nusselt numbers are lower than those for continuum flow and that 
they decrease with increased mean free path. The case of slip flow between parallel 
plates at equal temperatures is included for comparison; the length of duct necessary to 
obtain a fully developed Nusselt number is, for a given mean free path, the shortest for 
this case. 


INTRODUCTION 

The interest in the fluid- mechanics and heat-transfer characteristics of slightly 
rarefied gases flowing in conduits has been increased in recent years by the advances in 
high-vacuum technology and by the practical realization of flight at very high altitudes. 

At normal gas pressures the molecular mean-free-path length is so small that it is pos- 
sible to consider the behavior of the gas in a conduit by means of well-established 
theories of continuum fluid mechanics and heat transfer, with the gas adjacent to a sur- 
face assuming the velocity and temperature of the surface. At reduced pressures, on the 
other hand, such as might occur in a vacuum application, the gas flow may take place at 
a density low enough to permit a velocity slip and a corresponding temperature jump at 
the conduit wall. This density regime is called the regime of slip flow, and it is this 



regime with which the present investigation is concerned. In this regime the molecular 
mean-free-path length is still small but not negligible. 

The problem of laminar continuum flow between parallel flat plates with different 
temperatures prescribed along each of the two walls is considered in references 1 to 4. 
Heat transfer for laminar slip flow in a parallel-plate channel with constant wall temper- 
atures is studied in reference 5. The results allow for both walls at the same tempera- 
ture or else a uniform temperature at one wall and insulation at the other. 

Specific consideration is given herein to the heat transfer with laminar, incompres- 
sible, fully developed slip flow between a pair of parallel plates, the temperatures of 
which are maintained constant but not equal. A flat-duct heat exchanger, for example, 
which heats or cools a gas passing through it may be in contact with a constant tempera- 
ture bath at each plate. Since the baths are, however, at different temperatures, they 
thereby maintain the duct walls at constant but unequal temperatures. Another example 
of channel flow with constant but unequal wall temperatures arises when different coolants 
undergoing a change of phase are so employed in adjacent flow passages that the tempera- 
ture of one plate of the gas flow channel differs from that of the other plate. Such ex- 
changer channels are customarily made with small wall thickness and of material with 
high thermal conductivity, so that a wall is very nearly at a uniform temperature. 

The present analysis requires the calculation of the odd eigenvalues and eigencon- 
stants for laminar slip flow. The complete solution is then obtained by combining these 
odd quantities with the even eigenvalues and eigenconstants that have been determined in 
reference 5 for slip flow of a rarefied gas between two parallel plates with both walls at 
the same temperature. 

Since the convective term in the energy equation involves the velocity distribution, 
the first step in the analysis is to specify the velocity over the channel cross section. 

The fully developed gas velocity distribution for laminar slip flow in a parallel -plate 
channel was investigated in reference 6, and the results are used in the present study. 
With the velocity distribution specified, the energy equation can be considered. The so- 
lution involves series expansions, and asymptotic analytical expressions are derived for 
the eigenvalues and eigenconstants as functions of the gas rarefaction parameters and 
entrance-temperature parameter. Numerical solutions of these quantities are also ob- 
tained, and the two methods of computation are compared. 


ANALYSIS 

A schematic diagram of the parallel -plate channel showing dimensional nomenclature 
and coordinates is shown in figure 1. The direction of the gas flow is from left to right. 
For x < 0 the plates and the gas are assumed to be isothermal at temperature t g , 
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Figure 1. - Physical model and coordinate system for parallel-plate channel. 


whereas for x > 0 the upper plate 
temperature is step-changed to the 
constant value t . while the lower 

W, 1 

plate temperature is step-changed to 
the constant value t «. (All symbols 

W y £ 

are defined in appendix A. ) Each plate surface is characterized by a pair of properties, 
the thermal accommodation coefficient a and the reflection coefficient c r (ref. 7). In 
general, a and a may be different for each surface; for the present study, however, it 
will be assumed that a^ = ag s a and that cr^ = s a • Throughout the analysis, the gas 
flow is assumed hydrodynamically fully developed, viscous dissipation and axial conduc- 
tion are neglected compared with the conduction in the transverse y-direction, fluid prop- 
erties are assumed constant, and interactions between the velocity and temperature 
fields, such as thermal creep, are neglected. It is desired to determine the wall heat- 
flux distribution and the variation in the heat-transfer coefficient at each plate along the 
entire length of the channel. 


Energy Equation 

The starting point of the analysis is the differential equation for convective heat 
transfer in the parallel -plate channel flow and is 


pc u — = K 
p dx 



(i) 


The dimensionless velocity distribution u/ u, where u is the average velocity, is ob- 
tained from reference 6 and is given by 


u(v) = _ 3 (1 - 77 2 + 4a) 

u 2 1 + 6a 


(2a) 


S - f(l) 

u 


6a 

1 + 6a 


(2b) 


where a s £ u /2L. The slip coefficient | u is given by the expression (ref. 7) 

2 — 0" y) 


^u = 


(3) 


3 



where c r is the reflection coefficient and H is the mean free path: 


j t = 


12 p 


( 4 ) 


Equation (1), with u given by equation (2a), is to be solved subject to appropriate 
boundary conditions. Two conditions are obtained upon consideration of the effect of gas 
rarefaction on the thermal boundary condition at the two walls, which permit a jump be- 
tween the wall temperature t w and the adjacent gas temperature t^ (ref. 7): 





at r) = 1, x > 0 



at rj = - 1, x > 0 


(5a) 


(5b) 


where represents a temperature-jump coefficient related to other properties of the 
system by 




2 - a 2A l 
a X + 1 Pr 


( 6 ) 


The third boundary condition is a specified gas entrance temperature: 


t = t e at x = 0, - 1 < r\ < 1 


(7) 


Equations (1), (5), and (7) may be expressed in terms of nondimensional quantities as 


f (rj) 


0T _ d _ T 

\ 2 


( 8 ) 


T = T 


g, 


= 1 - at t? = 1 


= 1 , ? > 0 


(9a) 


T s T 


S> 2 


i-2r( d -I 




at r 7 = -1, £>0 


(9b) 
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where 


T s T 0 at 5 = 0, -1 < t? < 1 


t - t 

T _ w, m 

t i - t 
w, 1 w, m 


T = *g, 1 " t w, m 

S’ 1 " t ~t 

w, 1 w, m 


T = *g, 2 ' V m 

S, 2 " t “t 

w, 1 w, m 


T s e w, m 
e 

' ' " ' w, m 


V, 1 t - 


(9c) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


and 

r = s t /2L 

If a solution for T that will apply over the entire length of the channel is to be ob- 
tained, it is convenient to break T into the sum of two parts. The first part is 
T d s (*d ” Sy m^Sv i “ l w m^’ the fu Hy developed solution, which applies far down the 
channel from the entrance. The second part is T* = (t - ^/(t ^ - t m ), which is an 

entrance -region solution that is added to T^ to obtain dimensionless distributions in the 
region near the entrance of the channel. The dimensionless temperature distributions 
throughout the channel are given by 


T = T d + T* (14) 

Fully Developed Solution 

Far from the entrance to the channel, the gas temperature is independent of that 
is, it is a function of v\ and the gas rarefaction parameters only, and equation (8) yields 


5 



(15) 



= 0 


The solution of equation (15) is = A 77 + B, where A and B are constants. This 
result indicates that the temperature is linear across the channel. The boundary condi- 
tions are obtained from equations (9) as 


/9TA 

Tj(i) = 1 - 2r — 5 

v Vi 


T d (-1)=- 




(16a) 


(16b) 


For the fully developed situation, the boundary condition at the entrance of the heated 
channel (£ = 0) need not be considered, since it is accounted for by the entrance region 
solution. When equations (16) are used, the fully developed solution becomes 


T d = 


1 + 2 r 


V 


(17) 


Entrance Region Solution 

If the solution in the thermal entrance region is to be determined, the function T* 
is needed. The function T* must satisfy the equation 

(18 ) 

with the boundary conditions 


T*(l) = -2r/— ^ at r? = 1, ? > 0 

^ Vi 

(19a) 

T*(-l) = 2]V aT *'\ at 77 - -1, ? > 0 

(19b) 

CD 

■C* 

II 

1 
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At £ = 0, the condition is 


T(0,n) = T e = T d ( 77 ) + T*(0 ,rj) 


or, by rearranging, 


T*(0, t)) = T t) at ? = 0, -1 <n < 1 (19c) 

e 1 + 2r “ “ 

A solution to equation (18) may be obtained by the use of the separation-of-variable 
method, whereby T* is taken in the form of a product T* = X(£)H(?]), where X is a 
function of £ alone and H is a function of r\ alone. Substitution into the differential 
equation (18) shows that 


X = e-« 


( 20 ) 


2 

— + tfffajH = 0 
dry 2 


( 21 a) 


From equations (19), the boundary conditions on equation (21a) are 



( 21 b) 


Equations (21) comprise an eigenvalue problem of the Sturm- Liouville type. Solutions 
fitting the boundary conditions can be found only for discrete values of if>, that is, 
i $2 . • . ip n , termed eigenvalues. The solutions H^, ^2 * ' ' ®n 3X6 cor “ 
responding eigenfunctions. Hence, the solution for T* can be written as 


OO 



n=l 


( 22 ) 


The eigenfunctions H n are treated as functions of 17 , although, strictly speaking, they 
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are functions of both rj and the eigenvalues $ . 

The constants G n in equation (22) are evaluated to satisfy the condition at the en- 
trance to the heated channel (£ = 0). Evaluating equation (22) at £ = 0 and equating the 
result with equation (19c) gives 


oo 

Xj 

n=l 


G n H n (T 7) = T *(0’ ? ?) = T e 


l + 2r 


(23a) 


It follows immediately from the properties of the Sturm- Liouville system that 


G, 


n 



1 + 2T / 


??) f 0?) H n(??) dT 7 



f07)H„(T])d77 


*./ f(^)H n (l?)d7j 

f ' f(T?)Hj(7j)d7J 



1 

i + 2r 



V^(v)\(v) d V 

f(r?)H^(?7)d77 


(23b) 


Because of the symmetry of f ( 77 ) about the passage centerline, the eigenfunctions may be 
even or odd in 77 . If H n is an even function, the second term of equation (23b) disap- 
pears, while if H n is an odd function, the first term disappears. Hence, the constants 
G n divide into two classes a n and b n that are given by the results 


T e f f0?)Y n (T7)dTj 



n 




V f (v)Z n (v)dv 




(v) d V 


1 + 2r / 
«/o 


T? f (’?)Z n (r7)d7? 


J ^ f ( T 7) Z n( ?? )d V 


(24) 


(25) 
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where Y n (rj) and Z n (r]) are even and odd functions, respectively; that is, Y n (~v) = Y n (v) 
and Z n (-? 7 ) = -Z n (r/). The corresponding eigenvalues are denoted by 0 n and y n , and the 
solution for T* can be rewritten as 




-/3 £ 

^/n Y n^ e 

n=l 


♦I 

n=1 


b n Z n^ e 




(26) 


In reference 5 the result for the coefficients a n , when divided by the wall-temperature 
parameter T e , is shown to reduce to 



The series coefficients a„ are thus 

n 


a n = 


9Y + 2r 0 2 Y 



WP J v=1 .^ n 


(27) 


(28) 


This result implies that once the eigenfunctions Y n (?7;j3 n ) are known, the coefficients a n 
of equation (28) may be evaluated for a given value of T , for example, T = 1. For 
other values of T , it is only necessary to multiply a n (for T g = 1) by the new value of 
T e to obtain the new coefficients. 

As shown in detail in appendix B, the integral appearing in the numerator of equa- 
tion (25) may be written as 



T? f (’?)Z n (??)dr] = - 


1 + 2r 

^n 



(29a) 


whereas the denominator of equation (25) becomes 


9 


(29b) 



Hv) z liv) d v 



Then, the series coefficients b n are 


b = 2 — (30) 

\9y dridy/ * 

' r ‘i r/ T]=i ;y=r n 

Now that and T* are known, they can be superposed as in equation (14) to obtain 
the solution that applies over the entire length of the channel, which is 


T = 


i + 2r 


E -jS £ -y £ 

a n Y n^ e n + ) b n Z n ( ^ )e “ 
n=l n=l 


(31) 


Wall Heat Fluxes 


When the wall temperatures are specified, the wall heat-flux variations along the 
channel length required to maintain the wall temperatures constant are of practical in- 
terest. If a sign convention is adopted to give positive heat-transfer rates at both walls 
for £ — °°, then the heat-transfer rate from the upper wall to the gas is given by 


Vi 




and from the gas to the lower wall by 





(32) 


( 33 ) 


The heat-transfer rates q w ^ and q^ ^ can be f° un< ! by differentiating equa- 
tion (31) with respect to r\ and evaluating tiie result at t] = 1 and at 77 = - 1 to obtain 
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(34) 


OO oo 

q ,L , V' /dY\ ~P n Z /dZ \ -v 

w -’i = j— + > a ( — 2] e n +> b-JM e n 

*(V rVJ 1 + 2r '*> V ^ Vi 

n=l n=l 


n=l 

' ' ' Tj=l 

n=l 

OO 


OO 


■ V* e'^ + y b &) a'”- 

4w, 1 ‘ ‘w, m) 1 + 2r Z J n V*j/„ =1 Zj-V*!/,,.! 

n=l n=l 


(35) 


since 


f dY n\ _ / dY n\ 


‘ dr ?V! ^ d7 ? V-i 


and 


f dz n\ / di V 




Equations (34) and (35) indicate that the heat-transfer rates q w ^ and q w g become 
uniform and equal along the length of the duct after the development length, that is, after 
the influence of the eigenfunctions dies away. 


Nusselt Numbers 

In continuum heat-transfer theory, it is customary to represent the heat-transfer 
results in terms of a heat-transfer coefficient h = q w /(t w - t^) and a Nusselt number 
Nu h hD T A, where D,p is the thermal diameter and defined as 

j-j _ ^ Cross-sectio nal area 
T Heated perimeter 

For the flat duct, D T = 4L. The extension of these concepts to low-density flows is quite 
natural and appropriate. 

The Nusselt number for the upper wall may be written 
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(36) 


Nu, = 


_ V 1 4L 


V, 1 " *b « 


where t b , the mixed- mean temperature of the gas, is given by 

ut dr] 


t b _ 


/ 1 /•! 

rl 2 *^1 

udjj 


f(rj)t drj 


(37) 


since 


£' 


d v = 2d 


The Nusselt number Nu^ may be expressed alternatively as 


NUj = 


4 /ST\ 


where 


T = w, m = 1 r 

b t , - t „ 2 J-\ 


dr] 


w, 1 w, m 


(38) 


(39) 


The Nusselt number at the lower wall may be written in the same manner as 


Nu 2 4L 


(40) 


and hence 


Nu 2 = 


W 2 K 


4 / 9T\ 

1 + T b^V-i 


(41) 
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Putting equations (38) and (41) into more useful forms necessitates knowing the dimen- 
sionless bulk temperature T^, given by equation (39). Introducing the temperature 
distribution (eq. (31)) and the velocity distribution (eq. (2a)) into equation (39) yields 


T bijv v / mY » d ”iZ v ' Vn? / 

n=l J -\ n=l J -\ 


f(T?)Z n 


(42) 


The first integral on the right side of equation (42) can be evaluated by using equa- 
tion (21a) for an even eigenfunction with the result 



since ( < ^ Y n / < ^'>?) r7= o = 0 by symmetry, while the second integral in equation (42) reduces to 
zero, since integrals of odd functions from 77 = 1 to 77 = - 1 are zero. With these re- 
sults, equation (42) becomes 



(43) 


Final expressions for the Nusselt numbers may be written, with the use of equations (32) 
to (35), (38), (41), and (43), as 



4 


(44) 
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Equations (44) and (45) give the Nusselt number variation for each wall along the entire 
length of the duct. The fully developed value for the Nusselt numbers is obtained when 
£ — °° and is the same for both walls: 


Nu 


l,d 


= Nu 


2,d 


4 

i + 2r 


(46) 


Hence, it is necessary for the influence of all eigenfunctions (odd and even) to die away 
before this condition is reached. The fully developed condition also yields the linear 
temperature profile given by equation (17): 


T d = 


1 


i + 2r 


V 


(17) 


which remains unchanged at any value of £ from this point on. In the corresponding case 
with equal wall temperatures, the Nusselt number becomes fully developed when the tem- 
perature profiles at different distances from the entrance become similar, that is, when 
the value of £ is sufficiently large for the influence of the second and higher even eigen- 
functions to be negligible (ref. 5). This similarity occurs at distances which are much 
shorter than those in the unsymmetrical case. In addition, the fully developed condition 
in the case of unsymmetrical wall temperatures corresponds to the condition where there 
is no net heat transfer to the fluid, whereas in the symmetrical case there is a net heat 
transfer to or from the fluid. These factors result in (1) a shorter thermal entry length 
and (2) a higher value of Nu^ for the symmetrical as compared with the unsymmetrical 
wall temperature case. 

Equation (46) indicates that the fully developed Nusselt numbers Nu^ ^ and NUg ^ 
are dependent only upon the temperature -jump parameter T when the plates are main- 
tained at unequal temperatures. This result differs from that obtained when the wall tem- 
peratures are equal (ref. 5). In the latter case, the fully developed Nusselt number was 
found to depend upon the first eigenvalue given in reference 5, which, in turn, depended 
upon the two rarefaction parameters a and r. 
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Figure 2. - Fully developed Nusselt number ratio for slip 
flow in channel. Specular reflection coefficient, 1 ; 
ratio of specific heats, 1.4; Prandtl number, 0.73. 


Fully developed Nusselt numbers obtained 
from equation (46) are plotted in figure 2 as a 
function of p yR~t/2pL in the form of a ratio 
(Nu/NUq)^, where Nuq represents the con- 
tinuum value 4. 00. Also shown are fully de- 
veloped Nusselt numbers for the case of equal 
wall temperatures (ref. 5), where now Nuq 
= 7. 53. The effect of gas rarefaction is to de- 
crease the value of the Nusselt number below 
its continuum value. The effect on the Nusselt 
number of the specified value of accommoda- 
tion coefficient is apparent. The effects of gas 
rarefaction on the fully developed Nusselt num- 
ber are more pronounced with equal wall tem- 
peratures than in the corresponding case with 
unequal wall temperatures. 


Axial Variation of Gas Temperature Adjacent to Wall 

The axial variations of the gas temperatures adjacent to the wall are perhaps also of 
engineering interest. These quantities are obtained readily from equation (31) by setting 
r\ = 1 and t] = -1 therein. The gas temperature adjacent to the upper wall is 


oo op 

T<? ' 1)= T72? t Z a " Y '' <1)e " ,3n? + Z“ Z " 

n=l n=l 


-y n ? 
(l)e n 


while that adjacent to the lower wall is 


QO oQ 

T(? ’ _1> = ' I72T 

n=l n=l 


According to equation (21b), however, 


Y n (D 


= -2r 


(0Y„ 


ni 




(47a) 


(47b) 
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Z n (i) = -2r 



1 


Then, equations (47a) and (47b) become, respectively. 


T(S, 1) = — — 

i + 2r 



(48a) 


T(S,-1) = - 




(48b) 


In the absence of a temperature jump (r = 0), it is apparent from equation (48a) that 
t i = constant = t 1 , while from equation (48b) it is seen that for r = 0, t „ = t 0 . 

i z Wj ^ 

With a temperature-jump effect, however, t / t 1 and t 0 / t 

The foregoing expressions indicate that the solution requires the computation of jS^, 
y n , a n > t> n , an d (3Z n /3rj)^_j. Attention is now directed to the Sturm- 

Liouville eigenvalue problem (eqs. (21)). 


Transverse Distribution Functions Y(t)),Z(i7) 


The even function Y (rj) is the even solution of equations (21) and the customary nor- 
malization convention Y(0) = 1. Asymptotic expressions for large values of n for the 
even eigenvalues /3^ and constants a n (3Y n /3rj)^_j are given in reference 5, and 

the results are presented herein to make the analysis more complete: 


(4a + r) 


6 tan 6_ = 
n n 


V4« 


+ (1 + 4o!)sin - ^ 


4r(4a) 


3/2 


V l + 4 a )} 


(49) 



4ofT 


r(4a + r)(c + 1 + — 

V c / 


(50) 
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where 




I . t- . H . U . O A. 

Ratio of slip to average velocity, u s /U 

Figure 3. - Value of definite integral for any ratio of slip 
to average velocity. 


The first six roots of equation (49) are given 

in reference 8 for a number of values of C. 

The values of 1^ for any given ratio of slip 

velocity to average velocity u /u are shown 

s 

in figure 3. 

The odd eigenvalues y n and eigencon- 
stants b n and (dZ n /dj])^_^ are determined 
from the odd solution of equations (21) and 
the condition Z(0) = 0. This latter condition 
is obtained by setting 77 = 0 in the condition 
Z(?j) = -Z(-?7). If the methods presented in 
reference 5 to determine the asymptotic ex- 
pressions for large values of n for the case 


of equal wall temperatures are applied, the asymptotic solution of equations (21) for the 


odd functions Z(tj) satisfying the condition Z(0) = 0 is 


Z(?7) = D(1 + 4a!) 1 ' /4 (l - tj 2 + 4a)' 1 / 4 sin(yy J) (53) 

«)- 5/4 L S m(VTJ) t 4 Vvi(l-»l 2 t4 a ) :i/2 (Vy J ) co 3(Vy J )] (54) 
8rj 2 J Vl + 6a 

where D is an arbitrary constant and 

V3/^frj 

dT ? = - ■ 

The values of the function Z and its slope dZ/d n at the upper plate are found by setting 
r\ = 1 in equations (53) and (54), respectively; the results are 
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Z(l) = d (- - 4a ) 
\ 4a / 


sin e 


(56) 


1/4 



= — (1 + 4a) 1 / 4 



sin e 


2^ 1/2 (4a) 3/2 e cose 
J x (l + 6a) 1 / 2 


where 


e H V^ J i 


(57) 


(58) 


and 


J 


1 " 




(59) 


The eigenvalues y n are determined by the requirement Z(l) = -2r(dZ/drj)^_j. Then, 
with equations (56) and (57) combined in accordance with this relation, the eigenvalues 
are obtained as roots of the characteristic equation 


e n cote n = 


(4a + T) Via + (1 + 4a)sin _1 — 1 

L Vl + 4a _ 


4r(4a) 


3/2 


= -C (60) 

where e = J y The first six roots of equation (60) are given in reference 8 for a 
number of values of C. The eigenfunctions corresponding to these eigenvalues will be 
denoted by From equations (49) and (60) the eigenvalues j3 and y (and hence the 

eigenfunctions Y n and Z R ) can be seen to depend on the two rarefaction parameters 
£ U /2L and £ t /2L. 

If the series coefficients b fl (eq. (30)) are to be obtained, it is necessary to evaluate 
the terms (SZ/3y)^_ and (3 2 Z/3 rjdy) . These expressions are determined 

by differentiating equations (56) and (57), respectively, with respect to y. The results 
may be written as 
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(61a) 


(-) 

' 0r t?=i ;y=i 


D(1 + 4a) 1 / 4 C sin e 


n 


n 


2y n (4a) 


1/4 


9 2 Z 


D(1 + 4a) 


1/4 


C + 


(4a +r)(c +e 2 ) 




4y n (4a) 


cr 

5/4 


sin e 


n 


(61b) 


Thus, the coefficients b n as evaluated from equation (30) are given by 


b n = - 


8a(4a) 

w 


1/4 


D(4a + r)(l + 4a) 1 ' ’yC + 1 + J sin e n 


(62) 


The slope dZ/dr] evaluated at the wall, given by equation (57), can be alternately ex- 
pressed as 


'dZ\ D(1 + 4a) 1 / 4 sin e n 

v d? ^ rpl 2T(4a) 1 / 4 


(63) 


It is convenient to define a new coefficient B n , given by the product of b n and 
(dZ n /drj) ??=1 , as 



4a 


T(4a + r)lC + 1 + — 

\ c / 


(64) 


Equations (49), (50), (60), and(64)can be used to calculate the eigenvalues /3 n and y n and 
the coefficients A n and B n for arbitrarily chosen values of £ U /2L or u g /u and £ t /2L. 

The first four values of calculated from the asymptotic expressions together 

with the corresponding values of A n are listed in table I as analytical values for several 
values of the parameters u g /u and ^/2L. The first four values oi yfy n calculated 
from the asymptotic expressions together with the corresponding values of B n are 
listed in table n as analytical values for these same values of u g /u and |^/2L. The 
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TABLE I. - EVEN EIGENVALUES AND COEFFICIENTS FOR LAMINAR 


SLIP FLOW IN PARALLEL-PLATE CHANNEL WITH 


UNEQUAL WALL TEMPERATURES (T = 1) 




Ratio of slip to average velocity, u g / 

u 



1/3 

i 3/5 


Temperature- jump coefficient, l t /2L 


0. 1333 

0. 5333 

0.4 

1.6 





Solution 





Analyt- 

Numer- 

Analyt- 

Numer- 

Analyt- 

Numer- 

Analyt- 

Numer- 

I 

ical 

ical 

ical 

ical 

ical 

ical 

ical 

ical 

Eigenvalue 

#1 

1.431 

1. 168 

1. 289 

0.8128 

1. Ill 

0. 9087 

0. 8750 

0. 5273 

V^2 

4.321 

4.083 

4.010 

3. 630 

3.705 

3.602 

3.451 

3.337 

V^3 

7.265 

7.085 

6.906 

6. 682 

6. 645 

6. 586 

6.475 

6.417 

y/h 

10.28 

10. 14 

9.945 

9.794 

9. 700 

9.661 

9. 581 

9. 541 





Coefficient 




- A i 

0. 5245 

1. 292 

0. 1291 

0.6532 

0.4705 

0.8105 

0. 0869 

0. 2775 

” A 2 

.4455 

. 6554 

.0808 

. 1275 

. 1925 

.2077 

.0180 

.0196 

~ A 3 

.3380 

.4078 

.0441 

.0504 

.0789 

.0797 

. 00578 

. 00579 

~ A 4 

.2500 

. 2716 

.0255 

.0265 

.0399 

.0401 

. 00270 

. 00270 


even eigenvalues and coefficients for continuum flow (£ U /2L = £ t /2L = 0) in a channel with 
unequal wall temperatures are given in table HI and were obtained from reference 9, 
while the odd eigenvalues and coefficients are listed in table IV and were obtained from 
reference 4. 

The level of accuracy for the slip-flow eigenvalues and eigenconstants was checked 
by solving numerically the Sturm- Liouville equations (21) by means of the Runge-Kutta 
method on an IBM 7094 digital computer. Reference 5 presents numerical values for the 
even quantities K and A n calculated in this manner, and the results are given in 
table I. Equations (21) were solved numerically for the odd quantities in the course of 
the present investigation. The forward integration was started by using the condition 
Z n (0) = 0 and by arbitrarily letting (dZ n /d? 7 )^_Q = 1. The eigenvalues were found by 
trial and error until the boundary condition Z n (l) = - 2 T(dZ n /d 7 y)^_j was satisfied at 
77 = 1. The first four eigenvalues and coefficients B n calculated in this manner are 

given in table II. There is good agreement between the relevant quantities computed from 
the analytical asymptotic expressions and the numerical solutions even for values of n as 
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TABLE H. - ODD EIGENVALUES AND COEFFICIENTS FOR LAMINAR 


SLIP FLOW IN PARALLEL-PLATE CHANNEL WITH 


UNEQUAL WALL TEMPERATURES 


Ratio of slip to average velocity, u g / u 


1/3 


3/5 


Temperature-jump coefficient, 4^/2L 


0.1333 0.5333 0.4 1.6 


Solution 


Analyt- 

Numer- 

Analyt- 

Numer- 

Analyt- 

Numer- 

Analyt- 

Numer- 

ical 

ical 

ical 

ical 

ical 

ical 

ical 

ical 


Eigenvalue 


Vn 

2.871 

2. 610 

2.615 

2. 163 

2. 390 

2. 194 

2.051 

1. 849 

Vr 2 

5.780 

5. 574 

5.430 

5. 144 

5. 150 

5. 076 

4.950 

4.866 


8.755 

8. 603 

8.405 

8. 232 

8. 160 

8. 114 

8.025 

7.974 

VF 4 

11.79 

11. 68 

11.49 

11. 36 

11. 23 

11. 21 

11. 12 

11. 10 


Coefficient 


0. 4925 

0. 8717 

0. 1061 

0. 2503 

0. 3155 

0. 3868 

0.0398 

0. 0526 

.3925 

. 5099 

.0598 

.0761 

. 1199 

. 1231 

.00954 

. 00967 

.2921 

. 3290 

.0348 

.0354 

.0549 

.0547 

.00381 

. 00377 

.2139 

. 2244 

.0206 

.0201 

.0304 

.0298 

. 00201 

.00195 


TABLE m. - EVEN EIGENVALUES 
AND COEFFICIENTS FOR LAMINAR 
CONTINUUM FLOW IN PARALLEL- 


PLATE CHANNEL WITH UNEQUAL 
WALL TEMPERATURES (T 0 = 1) 
[Data from ref. 9.] 


Eigenvalue 


1. 372 

#2 

4. 625 


7.890 

#1 

11. 16 

Coefficient 

- A i 

1.708 

" A 2 

1. 138 

- A 3 

.951 

- A 4 

.848 


TABLE IV. - ODD EIGENVALUES 
AND COEFFICIENTS FOR LAMINAR 
CONTINUUM FLOW IN PARALLEL- 
PLATE CHANNEL WITH UNEQUAL 


WALL TEMPERATURES 
[Data from ref. 4.] 


Eigenvalue 

Vn 

2.992 


6, 250 

yjyi 

9.495 

yfyi 

12.78 

Coefficient 

B i 

1. 319 

b 2 

1.029 

b 3 

.894 

B 4 

.751 
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low as 3. The extent of agreement depends on the parameters u g /u and £^/2L with the 
best agreement at low values of 2L and higher values of u g /u. For those having 
access to digital or analog computing machines, the analytical expressions may provide 
helpful check solutions and for those who do not, the expressions presented herein em- 
ploying a rapidly converging series solution suitable for hand calculations and yielding 
values that are sufficiently accurate may be a useful tool. Furthermore, it is evident 
from the comparisons that the analytical expressions have definite application to deter- 
mining the higher eigenvalues for laminar slip flow heat transfer in channels. 


RESULTS AND DISCUSSION 

With the numerical information in tables I to IV, the variation along the length of the 
channel of the wall heat fluxes, Nusselt numbers, and temperature of the gas adjacent to 
each wall can be evaluated from equations (34) and (35), (44) and (45), and (48a) and 
(48b), respectively. These quantities also depend on the magnitude of the entrance tem- 
perature in relation to the values of the wall temperatures, or on the wall-temperature 
parameter T . 

The particular values of the wall-temperature parameter T g chosen for the compu- 
tations correspond to the following cases: (a) T g = +2, in which the wall temperatures 
are unequal but both are maintained at higher or lower temperatures than the tempera- 
ture of the gas entering the channel; (b) T g = Tl, in which the lower wall or upper wall, 
respectively, is at the temperature of the gas entering the channel; and (c) T g = 0, in 
which the arithmetic average of the wall temperatures is equal to the temperature of the 
gas entering the channel (i. e. , (t ^ + t 2 )/^ = t e ). 

The last case (T = 0) is a special case and warrants some additional discussion. 

For this special case, only odd eigenfunctions are involved in the various quantities of 
engineering interest (since a n = 0 from eq. (24)), and hence the wall heat-flux variations 
(eqs. (34) and (35)), bulk temperature (eq. (43)), Nusselt numbers (eqs. (44) and (45)), 
and temperatures of the gas adjacent to the walls (eqs. (48)) reduce to, respectively, 




*w, 2 J 




n=l 


= 0 
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Dimensionless wall heat flux, q w> jL/idt^ j ^ m ) 


i 



(a) Upper wall. Entrance temperature param- (b) Lower wall. Entrance temperature parameter, 0. 

eter, 0. 



Dimensionless axial distance, (x/2L)/RePr 

(e) Upper wall. Entrance temperature param- (f) Lower wall. Entrance temperature parameter, 2. 

eter, 2. 

Figure 4. - Heat-flux variation in thermal entrance region. Specular reflection coefficient, 1; ratio of specific 
heats, 1.4; Prandtl number, 0.73. 
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4 


Nu^ = Nug = 


^l +4 y 

l + 2r / y 
n=l 


V 


-y £ 
'n* 


and 


T(l) = -T(-l) = 


1 


i + 2r 


2r 


QQ 

Yb, 


-y £ 


L-j 

n=l 


The wall heat-flux variations required to maintain the wall temperatures constant 
were evaluated from equations (34) and (35) with the use of numerical values listed in 
tables I to IV and are plotted in figure 4 for several values of the rarefaction param- 
eter <p = n j/Rgt/2pL and the accommodation coefficient a and for the values of wall- 
temperature parameter T g of 0, 1, and 2. The results for T g = -1 and -2 are ob- 
tained from the figures for T g = 1 and 2, respectively, by interchanging subscripts 1 
and 2 on the dimensionless wall heat fluxes (e. g. , q w /c(t w j m ) f° r T g = -1 
is equal to q n L//c(t - t ) for T = 1). For almost all cases considered, an in- 
crease in gas rarefaction produces a decrease in the heat-flux requirement, at a given 
axial position, over that for continuum flow, for the range of the dimensionless axial dis- 
tance considered. The exception occurs for q w jLA(t w ^ - t w m ) with T g = 2. The 
accommodation coefficient a, which is generally near 1 for most surfaces and gases but 
which may take on much lower values for some gas-surface combinations, also has an 
important effect on the heat-flux variations along the channel length. Imperfect accommo- 
dation (a 4 1) may substantially decrease the wall heat-flux requirements because of the 
increase in effective thermal contact resistance associated with the increased tempera- 
ture jump. For the special case of T g = 0, the heat-flux variations are the same at the 
two walls, with heat flowing from one wall out through the other. 

The variation in Nusselt number along the channel length for each wall was evaluated 
from equations (44) and (45) for T g = 0, 1, and 2 with the use of the numerical eigen- 
values and eigenconstants from tables I to IV and is plotted in figure 5 for several values 
of the rarefaction parameter (p and accommodation coefficient a. The results for 
T = -1 and -2 are obtained from the figures for 
interchanging subscripts 1 and 2 on the Nusselt numbers. 

The effect of the gas rarefaction is, in general, to decrease the Nusselt number at 
each wall below its continuum value at every position along the channel. For most of the 
cases examined the Nusselt number decreases with increasing distance from the channel 
entrance; however, Nu^ for T g = 1 and 2 (and therefore Nug for T g = -1 and -2) 


T g = 1 and 2, respectively, by 
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Nusselt number along upper wall, Nuj 



(a) Upper wall. Entrance temperature parameter, 0. (b) Lower wall. Entrance temperature parameter, 0. 




Dimensionless axial distance, <x/2L)/RePr 


(e) Upper wall. Entrance temperature parameter, 2. (f) Lower wall. Entrance temperature parameter, 2. 

Figure 5. - Nusselt number variation in thermal entrance region. Specular reflection coefficient, 1; ratio of specific heats, 
1.4; Prandtl number, 0.73. 
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does not follow this behavior. In- 
finities and zeros occur in Nu, if 


1 < T e < 00 
_oo < T e < -1 


and in Nu 2 if 


Figure 6. - Variation of gas bulk temperature for flow in parallel-plate 
channel. Specular reflection coefficient, 1; ratio of specific heats, 
1.4; Prandtl number, 0.73. 


The Nusselt numbers 
become infinite when the bulk mean 
temperature of the fluid is equal to 
the appropriate wall temperature, 
that is, when = -1 or 1 (see 
eqs. (43) to (45)). The longitudinal 
variation of the dimensionless bulk 
temperature T^ (eq. (43)) was ex- 
amined and the results are plotted 
in figure 6. For T g = 1 and 2 the 
dimensionless bulk temperature T,. 
commences with the value of T g 
x = 0 and decreases with increasing 
distance from the entrance. For 


b 

at 


T = 0 

e 


the dimensionless bulk tem- 
perature remains constant (equal to zero) with increasing distance from the entrance. 

The effect of an increase in rarefaction is to increase T^ toward the entrance value T g 
at any given axial location. Imperfect accommodation increases this effect. The axial 
location, therefore, at which T^ = 1 or -1 moves downstream of the channel entrance 
as the gas mean free path is increased and/or the accommodation coefficient is decreased. 
The heat-transfer coefficient has little meaning for the situation in which the bulk mean 
temperature of the fluid equals the appropriate wall temperature. A thermal entrance 
length can be arbitrarily defined as the length of channel wall required for the Nusselt 
number to be within 5 percent of the fully developed value. The Nusselt number varia- 
tions Nu^/Nuj ^ and Nug/Nug ^ in the entrance region for unequal wall temperatures 
were evaluated as functions of the rarefaction parameter cp and accommodation coeffi- 
cient a and are presented in figure 7 for several values of the entrance temperature 
parameter T . A horizontal dashed line corresponding to an ordinate of 1. 05 or an or- 
dinate of 0. 95 (depending on whether the Nusselt number ratio approaches 1 from above 
or below) is shown in figure 7. The case of laminar slip flow between parallel plates at 
equal temperatures is investigated in reference 5 and the results are included (fig. 8) for 
comparison. In interpreting these figures, it is important to note that the Reynolds num- 
ber appears in the abscissa. For a given flow velocity and plate spacing, the Reynolds 
number is decreased with increased gas rarefaction. For continuum flow in a channel 
with unequal wall temperatures, the thermal entrance lengths are therefore very long, 
while the thermal entrance length at either wall is shortened with gas rarefaction. The 
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Nusselt number ratio along upper wall, Nu^/Nuj ^ 



(a) Upper wall. Entrance temperature parameter, 0. 



(c) Upper wall. Entrance temperature parameter, 1. 



(b) Lower wall. Entrance temperature parameter, 0. 

•o !- 6 r i \i i i i I l I I I “I r 


<\f 



I I I 1 l l l l I l i i I 

(d) Lower wall. Entrance temperature parameter, 1. 



Dimensionless axial distance, (x/2L)/RePr 


(e) Upper wall. Entrance temperature parameter, 2. (f) Lower wall. Entrance temperature parameter, 2. 

Figure 7. - Nusselt number ratio in thermal entrance region. Specular reflection coefficient, 1; ratio of specific heats, 1.4; Prandtl 
number, 0.73. 
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Dimensionless axial distance, (x/2L)/RePr 

(a) Upper wall. <b) Lower wall. 

Figure 8. - Nusselt number ratio in thermal entrance region. Specular reflection coefficient, 1; ratio of specific heats, 1. 4; 

Prandtl number, 0.73; equal wall temperatures. 

increase in temperature jump with a decrease in accommodation coefficient also has an 

important effect on the thermal entrance length. For a given mean free path, the length 

required for the Nusselt number ratio to be within 5 percent of its fully developed value is 

least for the channel with the walls at the same temperature (fig. 8). 

The variation along the length of the channel of the temperature of the gas adjacent to 

the upper wall and to the lower wall has been evaluated from equations (48a) and (48b), 

respectively, and plots are given for T g = 0, 1, and 2 in figure 9. The results for 

T = -1 and -2 are obtained from the figures for T =1 and 2, respectively, by in- 

terchanging subscripts 1 and 2 on the gas temperatures t^ and then multiplying the 

results by -1 (e. g. , T(£, 1) for T g = 1 is equal to -T(£, -1) for T g = -1). 

In the absence of rarefaction effects, the temperature of the gas adjacent to the walls 

changes discontinuously from the value t g at the channel entrance (x < 0) to the value 

t . or t _ at the appropriate wall and remains constant for x > 0. In other words, 
w, 1 w, 2 — 

the difference between the surface temperature t 7 . (or t 9 ) and the contiguous gas 

W y X lAf y & 

temperature t g ^ (or t g 2 ) is zer0 along the entire length of the channels. This result 
means that, for’continuum flow, T(£, 1) = 1 and T(?,-l) = -1 for x> 0 for all T g . 

With rarefaction, the temperature differences t - t w ^ and t^ 2 " *w 2 have non_ 
zero values. In general, the effect of an increase in rarefaction is to decrease T(?, 1) 
and increase T(£, -1) toward the entrance value T g at any given axial location. The in - 
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crease in the difference between the local wall temperature and the contiguous gas tem- 
perature with an increase in accommodation coefficient is evident. The solutions for 
T(£, 1) and T(£, -1) approach their respective limiting values at x-«° asymptotically 
for the slip flow of a gas, while for the continuum flow, there is, as mentioned pre- 
viously, a step change in T(£, 1) and T(£, -1) at x = 0. 


CONCLUSIONS 

An analytical study was conducted for the laminar forced- convection heat transfer to 
a slightly rarefied gas flowing between two long parallel surfaces. These surfaces were 
taken to have constant but unequal temperatures. The wall heat fluxes and Nusselt num- 
bers in both the entrance and the fully developed regions were obtained as functions of 
the gas mean free path for various values of the wall temperature constant T . The 
analysis provided analytical expressions for the even and the odd eigenvalues and eigen- 
constants which show good agreement with numerical results, especially for the higher 
eigenvalues and eigenconstants. Some of the characteristics of laminar slip flow heat 
transfer can be summarized as follows: 

1. The wall heat-flux variation required to maintain the wall temperatures constant 
is, in general, decreased with an increase in gas rarefaction and/or a decrease in the 
value of the accommodation coefficient. The heat flux at either wall approaches asymp- 
totically a constant value that depends solely on the temperature-jump coefficient. 

2. The Nusselt numbers do not become constant until a linear temperature gradient 
is established in the gas, which requires, for continuum flow, fairly long entrance 
lengths. The Nusselt number, in general, is decreased below its continuum value because 
of the effect of gas rarefaction and imperfect thermal accommodation. 

3. The fully developed Nusselt number is the same for either wall and is independent 
of the magnitude of the entrance temperature in relation to the values of the wall tempera- 
tures. It does, however, depend on the temperature- jump coefficient. The effects of gas 
rarefaction on the fully developed Nusselt number are more pronounced in a channel with 
equal wall temperatures than in a duct with unequal wall temperatures. 

4. The thermal entrance length at either wall is shortened from its continuum flow 
value with gas rarefaction. For a given mean free path, the thermal entrance length is 
least for the channel with the walls at the same temperature. 

Finally, other rarefaction effects, such as thermal creep velocity, would undoubtedly 
moderate the results given herein. Such effects were not considered in this analysis. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 14, 1965. 
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APPENDIX A 


SYMBOLS 


A n 

coefficient defined by eq. (50) 

Nu 

Nusselt number, h D ^/k 

a 

accommodation coefficient 

Pr 

Prandtl number, (xc^/k 

a n 

even coefficient in series ex- 

P 

gas pressure 


pansion 

^w 

rate of heat transfer per unit 

a n 

even coefficient divided by 

area from wall to gas 


T g , i. e. , a n /T e 

Re 

Reynolds number, 2puL/p 

B n 

coefficient defined by eq. (64) 

R_ 

g 

T 

gas constant 

b n 

odd coefficient in series ex- 
pansion 

dimensionless temperature, 
- *w, 1 - ‘w, m> 

C 

constant defined by eqs. (49) 
and (60) 

T* 

dimensionless entrance region 
temperature, (t - t^)/ 

C P 

specific heat of gas 


^w, 1 “ *w, 

D 

arbitrary constant 

t 

gas temperature 

d t 

thermal diameter for channel, 
4L 


temperature of gas adjacent to 
wall 

f(i?) 

dimensionless velocity, u(rj)/u 

t 

w, m 

temperature equal to arithmetic 

G 

coefficient in series expansion 


average of wall temperatures, 

n 

H(tj) 

h 

function of 77 

heat- transfer coefficient, 

u(?7) 

K, 1 * V 2>/ 2 

gas velocity 


~ V 

u 

average velocity 

h 

definite integral defined in 

x«) 

function of £ 


eq. (52) 

X 

axial coordinate 

J 

indefinite integral defined in 
eq. (55) 

Y 

even transverse distribution 
function 

J 1 

definite integral defined in 

Y n 

even eigenfunctions of eqs. ( 21 ) 


eq. (59), equal to 1 ^ 

y 

transverse coordinate 

L 

half distance between plates 

z 

odd transverse distribution 

SL 

mean free path 


function 
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Z n odd eigenfunctions of eqs. (21) 

Z v odd eigenfunctions of eq. (B 3 ) 

a dimensionless velocity- slip 

coefficient, l u /2L 

/3 separation constant 

/ 3 n even eigenvalues of eqs. (21) 

r dimensionless temperature-jump 
coefficient, £ t /2L 

Y separation constant 

y odd eigenvalues of eqs. (21) 

y v odd eigenvalues of eq. (B 3 ) 



£ dimensionless axial distance, 

4 (x/ 2 L)/RePr 

77 dimensionless transverse coordi- 

nate, y/L 

k gas thermal conductivity 

X ratio of specific heats 


M gas viscosity 

(jj. temperature- jump coefficient 

l u velocity-slip coefficient 

p gas density 

a specular reflection coefficient 

<p rarefaction parameter, p ^Rgt/2pL 

separation constant 
Subscripts: 

b gas bulk condition 

d fully developed region or condition 

e entrance, x = 0 

g gas adjacent to wall 

s slip 

w wall 

0 continuum flow conditions 

1 upper wall, y = L 

2 lower wall, y = -L 
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APPENDIX B 


EVALUATION OF SERIES COEFFICIENTS b n 

The odd coefficients of the series expansion (eq. (26)) axe determined as the quotient 
of two integrals (eq. (25)): 


n 


i + 2r •'o 


T?f(7?) Z n 07) d ?7 



f(r?) z n(i?) d i? 


(25) 


With equation (21a) applied to the odd eigenfunctions Z n (rj), the integral appearing in the 
numerator may be written as 



since Z n (0) = 0. According to equation (21b), however, 

z " (1) ' - 2r ©„ =1 

so that 



?7f(T?) z n (T?) d T7 = " 


1 + 2T 


y n 



1 


(Bl) 


(B2) 
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In order to evaluate the integral in the denominator of equation (25), let and Z n 
be the solutions associated with two distinct values of y: y^ and y^. Therefore, 


and 


d 2 Z 


&n‘ 


- + yf(ij)Z„ = 0 with Z (1) = -2r 


^dZ 


d7 ? Vi 


; z (o) = o 


d 2 Z /dZ \ 

5 + y n f(T ? )Z n = 0 with Z n (1) = - 2 r ( — ) J z n (°) = 0 

dr] 2 VdT? V 1 


(B3) 


(B4) 


Equation (B3) is multiplied by Z n and equation (B4) by Z p , and then equation (B4) is 
subtracted from (B3). The result, after transposing, is 


A l dZ dZ ) 

(y - y )f ( t?)Z Z = - — (Z — - - Z 

lr ” v " n d„ \ " d „ n d„ y 


If this equation is integrated between 0 and 1, the following is obtained: 


f l ( dZ n dZ „\ 

<v„-W 


(B5) 


where Z v (r],y v ) and Z n (rj,y n ) for v ^n are orthogonal functions with respect to the 
weight function f ( 77 ); that is, 



for y / n 


(B 6 ) 


This property of the eigenfunctions was used in obtaining equation (25). 

If y v = y n , the integral on the left side of equation (B5) becomes the integral of the 
square of the characteristic function: 
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The expression on the right assumes the intermediate form 0/0. Hence, it is necessary 
to apply L'Hospital’s rule and to differentiate numerator and denominator with respect to 
y v before setting y^ = y n> Carrying out this differentiation results in 



f(^)Zn(r})dr] = 




r 3z 9Z„ 
n n 

9 y n/ 


rj=l 



az„ 9Z 

n n 


9y. 


(B7) 


n, 


rj = 0 


Near rj ~ 0, the boundary condition Z n (0) = 0 requires that the odd solution to equa- 
tion (21a) be given by Z n (T] ~ 0) = F sin 0 )tj^, where F is an arbitrary constant. 

Hence, at rj = 0, 





= 0 


Thus, equation (B7) reduces to 



Uj])Z^(ri)dr} = 



9Z 


n 


dy i 


- Z 


n 


n 


9 2 Z. 


n 


dn dy 


ri=l 


(B8) 


Into this expression, there is now substituted the boundary condition 


z n (i) = -2r 



T]=l 


This substitution gives 
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f(r?)Z^O})d77 



+ 2r 


d 2 Z \ 

dr i r=r n 



Therefore using equations (25), (B2), and (B9) yields 


b 


n 


1 



+ 2T 


3 2 Z \ 

dri %=i;^ 


(B9) 


(BIO) 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere atjd space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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